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Алгоритм расчёта величины износа на поверхности зуба во время приработки 
цилиндрической косозубой передачи” 


А. Ю. Луконин 


При приработке передачи происходит износ зубьев. Описан алгоритм расчёта величины такого износа в ко- 
созубой цилиндрической передаче с точечным контактом. Работа алгоритма демонстрируется на примере 
расчёта косозубой цилиндрической передачи с круговым профилем зуба в торцевом сечении. Задача реша- 
ется с помощью программной системы конечноэлементного анализа АМ$У$. Описано применение разрабо- 
танного алгоритма для конечноэлементной модели зацепления пары зубьев зубчатой передачи. Проведён 
анализ изменения поверхностного распределения контактных давлений в зависимости от величины износа, 
полученного при приработке зубьев. Установлено следующее. Зубчатые передачи с круговым контуром зуба 
в торцевом сечении обладают новыми свойствами по сравнению с передачами, имеющими точечный или 
линейный контакт. Рассмотрена приработка зубьев в зацеплении с круговыми контурами. В этом случае при 
увеличении износа максимальные контактные давления уменьшаются и перемещаются со средней линии 
зуба в центр пятна контакта. 
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Введение. В процессе эксплуатации зубчатых передач происходит износ рабочих поверхностей 
зубьев. Он может по-разному влиять на несущую способность передачи. Иногда износ приводит к 
разрушению контактирующих поверхностей (так всегда происходит, например, в эвольвентном 
зацеплении). Но износ может способствовать и улучшению контактной несущей способности пе- 
редачи — например, в процессе приработки передач с зацеплением Новикова [1] (износ типа ис- 
тирания [2]). Аналогичный процесс происходит при приработке передач эллиптического типа [3] 
с круговым профилем зуба в торцевом сечении. 

При нагружении зацепления эллиптического типа в локальных зонах контакта на средней 
линии зуба возникают большие контактные давления. В процессе приработки поверхностей они 
вполне могут нивелироваться в силу высокой степени согласованности поверхностей контакта. 
Поэтому в данной работе исследуется процесс приработки поверхностей контакта эллиптического 
зацепления на начальном этапе эксплуатации зубчатой передачи. 

Алгоритм оценки величины износа на рабочей поверхности зуба. Рассмотрим процесс 
приработки передачи, которая происходит в результате истирания поверхностей контакта. Тогда 
можем предположить, что величина (г износа на площадке &а5 поверхности пропорциональна 
работе А сил трения на этом участке [4]: 

А=В,,-а!- 45, (1) 
где В, — константа, определяющая свойства материала. 


При вращении зубчатых колёс пятно контакта в зацеплении эллиптического типа движет- 
ся вдоль средней линии зуба. Поэтому работу сил трения на площадке 46$ за один цикл контакта 
можем записать: 


ь 
А= а5К, | ИВ (ЕЕ, (2) 
Е 


"Работа выполнена в рамках инициативной НИР. 
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где У, — скорость скольжения; К, — коэффициент трения; р(Ё) — контактное давление на 


площадке &@5; Ё — момент времени вхождения площадки &@$ в пятно контакта; Ё› — мо- 


мент выхода из пятна контакта. 
Скорость скольжения в пятне контакта постоянна: 


И =Г(@, +0, ), (3) 
Здесь г — радиус кривизны головки и основания зуба (рис. 1), ®,, х, — соответственно 


угловые скорости вращения колеса и шестерни. 
У 





Рис. 1. Величина износа 


Поэтому скорость скольжения можно вынести из-под знака интеграла и перейти к инте- 
грированию по статическому пятну контакта. Для этого сделаем замену переменных 


а! =У Е, (4) 


где У, — скорость пятна контакта вдоль средней линии зуба, которая является постоянной вели- 


ЧИНОЙ. 
Тогда выражение (2) для работы А за один цикл контакта можем записать в виде 
Ки &5 
_ `СТро ск 
И | р, (191, (5) 
где /[ — кривая, параллельная средней линии зуба и пересекающая центр площадки @&@5; 


р, (Г) — контактные давления на этой кривой. 


Из (1) и (5) для величины износа зуба в любой точке кривой [ получаем 


К, У 
т |2. (а. (6) 


мат п (Е 





В каждом торцевом контуре зуба построим полярную систему координат, как показано на 
рис. 1. Из (6) следует: для любой точки с полярной координатой а на контуре зуба зубчатого ко- 
леса, вращающегося с угловой скоростью и, получаем выражение для вычисления величины из- 
носа за время Ё 

ки 
а, (Е) =С(Ё)9., 9. = | Р. (14, С = у ФЕ, (7) 
р, 


мат п 
где /[, — кривая [, проходящая через точку с полярной координатой а в полярной системе коор- 
динат контура торцевого сечения (см. рис. 1). 
Из (7) получаем, что если обозначить 4 (Е =Д,, то для величины износа в точке с про- 


извольной координатой а можем записать: 
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ат. = ЗА. (8) 
9 


Оценка величин 9, проводилась с помощью системы конечноэлементного анализа АМ$У$. 


Конечноэлементная модель зуба для исследования процессов приработки его поверхности пред- 
ставлена на рис. 2. 
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Рис. 2. Конечноэлементная модель зуба 


Для оценки о выделим на поверхности контакта слой конечных элементов ВДОЛЬ ЛИНИИ 


[, шириной В (см. рис. 2). Тогда согласно (7) можем записать 
1 
[ер Р ь2.Р,5, Г. (9) 


где суммирование проводится по элементам выделенного слоя, $, — площадь го контактного 


поверхностного элемента, р, — контактное давление в Ем элементе. 


Выражения (8) и (9) описывают расчёт приработки зуба, в предположении, что контакт- 
ные давления не зависят от времени. Это справедливо при малых Ё{и аг.. Если это допущение не 
выполняется, то соотношения (8) и (9) из выражений для вычисления величины износа превра- 
щаются в систему нелинейных уравнений для вычисления величин д". (Е) и р, (а^. (Е). 

Согласно общей методологии исследования нелинейных систем в АМ5У$ для решения этой 


нелинейной системы уравнений использована пошаговая процедура расщепления нелинейных 


связей. 
[Шаг 0. 


Задаём начальную форму поверхности контакта с помощью узловых точек 7) ‚ Здесь /— 


номер узла в }-м торцевом сечении зубчатого колеса. При этом узлы равномерно (с шагом БЬ) рас- 
пределены вдоль контура каждого торцевого сечения на поверхности контакта, а множества уз- 


лов РО ри Ро) задают границы го слоя конечных элементов (рис. 2). В начальной конфигура- 


ции зуба в торцевых сечениях узлы лежат на дугах радиуса г (радиус головки зуба) [3]. 
Цикл по шагам А =1,..., М. На каждом шаге процедуры предполагается, что в точках ба- 


зового слоя г“) величина износа равна Д.. 
Шаг ^. 
1. Используя координаты узлов 7“) на контактной поверхности, построим конечноэле- 
К У 
ментную модель зацепления и рассчитаем поле р“) контактных давлений (/— номер контактного 
элемента на поверхности контакта). 


2. Вычисляем 9“ с помощью (9). 
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3. Используя (8), вычисляем величины приработки в узлах контактной поверхности: 





(К) 
(К) _ 9! 
АГ, ^ = Се) А 
0 
4. Находим новые координаты узлов, которые они получают в результате вычисления 
приработки зуба на А-м шаге процедуры: 
(К) _ (КУ (К) 5 (К-И 
Ио =Г, — аг, Пр у 


—>(К-—1 —(К-1 
где п, Е единичная внешняя нормаль к поверхности контакта в узле с координатами Г . 


Конец Шага (. 
Изменение контактных давлений в процессе приработки зубчатой пары. Оценка поля 
контактных давлений в эллиптическом зацеплении проводилась для зубчатой пары, нагруженной 
так, что окружная сила в зацеплении равна 8000 Н. При этом радиус головки зуба равен 3 мм. 
Контактные давления на рабочей поверхности исходной конфигурации зуба, рассчитанные в си- 
стеме АМ$У$, представлены на рис. 3. 





Рис. 3. Распределение контактных давлений на рабочей поверхности зуба 


Максимальные контактные давления в выбранном примере достигают 754 МПа в области 
точки А (см. рис. 3), а в области точки Суровень давлений не превышает 130 МПа. 

Для исследования процесса изменения контактных давлений в процессе приработки пере- 
дачи к данному примеру нагружения зубчатой пары был применён алгоритм расчёта износа, опи- 
санный выше. На каждом шаге итерации в качестве базового выбирался слой вдоль средней ли- 
нии зуба В — А — В (см. рис. 3). Приращение Д, величины износа в базовом слое устанавлива- 


лось равным 10`° м. В результате приработки зуба до величины износа в базовом слое, равной 
10? м, происходит изменение контура зуба, показанное на рис. 4. 





Рис. 4. Изменение контура зуба при приработке передачи 


Распределения контактных давлений при величине износа в базовом слое, равной 
5.10 © ми 10 > м, представлены на рис. 5. 
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а) 


6) 


Рис. 5. Изменение контактных давлений в процессе приработки передачи при величине износа в базовом слое, равной 
5.108м (аи 10° м (6) 


Точками А на рис. 5 обозначены зоны максимальных давлений. ИЗ рис. 3 и 5 следует, что 
при увеличении износа пятно контакта разбивается на два не связанных пятна, а область макси- 
мальных давлений перемещается в центр каждого пятна. При этом максимальные давления 
уменьшаются с 754 МПа при отсутствии износа до 501 МПа при величине износа в базовом слое, 
равной 10” м. 

Выводы. Проведённые исследования показывают, что контакт поверхностей в зацеплении эл- 
липтического типа обладает новыми свойствами, по сравнению с герцевским контактом. В зубча- 
тых парах с контактом Герца приработка ассоциируется лишь с процедурой шлифования контак- 
тирующих поверхностей. А в процессе приработки контакта в зацеплении эллиптического типа 
существенно меняется распределение контактных давлений по поверхности контакта. Кроме того, 
в этом случае снижается уровень максимальных давлений. Так, в рассмотренном примере при 


величине износа - 10? м максимальные контактные давления уменьшаются более чем в 1,5 раза 
по сравнению с давлениями в зацеплении с исходной геометрией. Данное свойство является 
следствием большей согласованности поверхностей контакта в эллиптическом зацеплении по 
сравнению с контактом Герца. 
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